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ABSTRACT 

Motivated by recent simulations of galaxy formation in which protogalaxies acquire their baryonic 
content through cold accretion, we study the gravitational fragmentation of cold streams flowing into 
a typical first galaxy. We use a one-zone hydrodynamical model to examine the thermal evolution of 
the gas flowing into a 10 8 M Q DM halo at redshift z = 10. The goal is to gain an understanding of the 
expected fragmentation mass scale and thus the characteristic mass of the first population of stars to 
form by shock fragmentation. Our model accurately describes the chemical and thermal evolution of 
the gas as we are specifically concerned with how the chemical abundances and initial conditions of 
the low density, metal enriched, cold accretion streams that pass an accretion shock alter the cooling 
properties and tendency to fragment in the post-shock gas. Cold accretion flows are not shock heated 
at the virial radius but instead flow along high-baryonic-density filaments of the cosmic web and 
penetrate deep into the host halo of the protogalaxy. In this physical regime, if molecular cooling is 
absent because of a strong LW background, we find there to be a sharp drop in the fragmentation mass 
at a metallicity of Z ~ 10 -4 Z Q . If, however, H 2 and HD molecules are present, they dominate the 
cooling at T < 10 4 K, and metallicitiy then has no effect on the fragmentation properties of the cold 
stream. For a solar abundance pattern of metallicity, O is the most effective metal coolant throughout 
the evolution, while for a pair instability supernova (PISN) metallicity yield, Si + is the most effective 
coolant. PISN abundance patterns also exhibit a slightly smaller critical metallicity. Dust grains are 
not included in our chemical model, but we argue their inclusion would not significantly alter the 
results. We also find that this physical scenario allows for the formation of stellar clusters and large, 
10 4 Mq bound fragments, possibly the precursors to globular clusters and supermassive black holes. 
Lastly, we conclude that the usual assumption of isobaricity for galactic shocks breaks down in gas of 
a sufficiently high metallicity, suggesting that metal cooling leads to thermal instabilities. 
Subject headings: cosmology: theory — galaxies: formation — galaxies: high-redshift — stars: for- 
mation hydrodynamics — intergalactic medium 



1. INTRODUCTION 



An important open question in cosmology is to iden- 
tify the conditions under which the first Population II 
(Pop II) stars could have formed. It has been argued 
that clusters of predominantly low-ma ss stars could have 
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dies 


the end of the so-called 



10 ( |Bromm fc Clarke 
Their emergence strad- 



the most uncharted time in the history of th e universe 
( |Barkana fc Loeb||2001| |Miralda-Escude||2003[ ). Between 
the release of the cosmic microwave background (CMB) 
at z ~ 1, 100 and t he highest redshift gamma ray burst 



(GRB) at z ~ 8.2 (Salvaterra et al. 2009 Tanvir et al 



20091, no direct observations have been made yet, thus 



rendering theoretical models our best tools for exploring 
this formative cosmic epoch (for indirect pro bes of this 
era, s eejGispert et al.|[2000l IChristlieb et al.| [2002{ |Cen] 



2003) |Frebel et al.||2005l |2009^ ~ 

The standard cosmological model (ACDM) predicts 
c first stars (Pop III) formed at redshifts z ~ 20 — 30 



in ~ 10 6 M 


Pt dark matter 'minihalos' (e.g., Couchman & 


Rees|1986 


Haiman et al.|1996 


Tegmark et al.|1997|) and 


this brought about an end to t 


le cosmic dark ages. The 



cooling properties of H 2 , the most effective primordial 
gas coolant at temperatures below 10 4 K, set a charac- 



teristic mass scale of 100 A fe for Pop III stars (Abel et al. 
2000 Bromm et al. 2002). The primordial initial mass 



function (IMF) would then have been top-heavy. This 
result has been confirmed by subsequent studies (e.g., 
l Omuka i & Pallaj |2003| IQ'She'a fc Norman|[2007l |Yoshida 
|et al.|2006[|2008[ ). Further investigations were directed at 
elucidating their proper ties and key implications, such as 
the chemical signature ( |Heger fc Woosley|2002| Jlwamoto 



et al. 2005), the detailed IMF ([Nakamura & Umemura 
20011) , the effect on reionization f Wyithe fc Loebl|2003[ 



I Sokasian et al. 20041, the formation rate ( jlrenti ez Sti- 



aveffil 2009) 



eFaT 2009 



and the tendenc y toward bmarity ([lurk 
Stacy et al. 2010). These massive Pop 

10 5 K 



sec 



stars had extremely hot photospher es (T e g 
Bromm et al.pOOla Schaerer|2002 ) , and very short life 



times of 
picture see 



Clark et al. 



3 Myr (for the case against this standard 
Tumlinson et~ai^|2004l |Silk fc Lang"er]|2006] 



2010p . 



It is still not known with any certainty when and how 
the universe shifted its star formation mode from the 
top-heavy IMF for P op III to the roughly S alp eter dis- 



tribution seen today. [Bromm et al. (2001b I suggested it 
is the chemical input from fop 111 stars that allowed gas 
in minihalos halos to cool more efficiently and fragment 
at a smaller mass scale. They found that when molecular 
cooling is excluded by a soft-UV photodissociating back- 
ground, gas enriched to Z cr it — 10~ 3 5 Z Q , the so-called 
critical metallicity, fragments into smaller mass clumps. 
The idea of metallicity being the main driver of the mode 
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transition was advanced by many other studies. ISantoro] 
& Shull (2006) found individual critical metallicities for 
C, Ke, Si, and O and suggested that for densities and 
temperatures conducive to Pop III star formation, these 
species contribute more to the c ooling rate than atomic 



a galaxy, leading to shocks at radii much smaller than 



the halo y irial radius (see, also, |Wise fc Abel 
ney l|2004| ). 
hat thi 



This differs qualitatively trom ho 



2007 Bin 



hydrogen or H 2 . ISmith et al. (12009) found that there 



exists a critical metalhcity window, due in part to the 
CMB temperature fl oor, which regulates the transition. 
In a s eries of papers ( Omukai||2000 Omukai et al. 112005" 



accretion 

in that the infalling gas passes through the virial shock 
unaffected, has a much higher density, and is deposited 
much closer to the center of the galaxy. Cold accretion 
is theoretically very appealing. It presents a solution to 



2008 1 dust was shown to have a significant impact on 
the Pop III /II transition. Its inclusion lowers the critical 
metalhcity derived from pure metal line cooling by sev- 
eral orders of magnitude. Dust grains act as a catalyst for 
molecule formation and radiate thermally, however these 
effects generally become important late in the pre-stellar 
gas contraction p hase when density is high. Schneider] 
fc Omukai ( 2009 1 argued that the interplay of dust and 
metaf cooling is the key driver, together with the CMB 
temperature floo r. A different view was deve loped in a 
series of studies ( Jappsen et al. |2007 2009a b), now tak- 
ing into account molecular cooling. According to these 
authors, metalhcity plays no role in the fragmentation of 
collapsing clouds while initial conditions, set by the de- 
tails of galaxy formation, the CMB temperature, turbu- 
lence, and rotation, dominate the Pop. III/II transition. 
This star formation mode transition is clearly a complex, 
gradual process, regulated in part by all of these factors. 

Most studies of the critical metalhcity to date have 
focused on ~ 10 6 M DM halos, or minihalos, as the 
primordial star formation environment. However, the 
majority of the first Pop II stars likely formed in deeper 
potential wells, possibly in halos that can coo l via atomic 
hydrogen cooling (e.g. 



the o bserved blue-red galaxy bimodality (Baldry et al. 
2004 1 . It provides a natural explanation for high red 



shift galaxies with large star formati on rates (SFRs) bu t 
no evidence of a major merger (e.g. JGenzel et al.|2006] >. 
Also, cold accretion naturally arises in many simulations 
of galaxy formation flKeres et al.|2005||Dekel fc Birnboim 
20061 1 Wise & Abel|2007| i Greif' et al.||2TO8l |Harford et al 



2008| |Regan fc Haehnelt 2009p , suggesting this type of 



accretion is natural to high-cr peak halos at high redshifts 
within ACDM models. 

DM halos with masses ~ 10 8 Mq represent the small- 
est structures that can be classified as 'first galaxies' (|Ki- 
" Read et all|2006l IWise fc Abe 



tayama fe Yoshida|[2005l [B 
2007| |Creif et al. 12008T |B 



this mass have virial temperatures T v j 



romm et al. 2009[ ). Halos of 



> 



10 4 K which 



allows efficient cooling via excitation of atomic hydrogen 
e.g., | Haiman et al.||1997[ |Tegmark et a l. [l997[ |Miralda - 

e & Rees 



Escui 



1998; Barkana & Loeb 



1999; 



rOh|[1999 



r 



This is the halo mass scale where we first expect to see 
sustained, possibly self-regulated, star formation, similar 
to the present-day case. Finally, in 10 8 Mq halos, gas ac- 
creting a long cosmic filaments will drive supersonic tur- 
bulence (IWise fe Abel 2007;|Greif et~aI1[2008l |Regan fe 



Haehneltp009[ IKlessen fe Henne bellc 2009), fun damen- 



tally changing the star formation dynam i cs (|Padoan fc 



Bromm et al. 2009). Whereas [Nordhmd 2002; Mac Low & Klessen 2004; McRee fc Os 



the quasi-hydrostatic, roughly spherical nature of bary- triker||2007p. In this work, as in many others, we define 



onic collapse inside of minihalos admits modeling of low 
metalhcity cooling in a rather straightforward fashion, 
the hydrodynamics inside the first galaxies is much more 
complex. Simulations with realistic cosmological initial 
conditions have recently demonstrated how cold accre- 
tion streams feed dense, turbulent gas to the center of 
the galaxy, wher e second generation star formation will 



take place (e.g., Greif et al. 2008 Wise fc Abel 2007 



20081. Since detailed simulations of the fragmentation 



properties in such a cloud with different levels of pre 
enrichment are still lacking, it is timely to carry out an 
exploratory survey of the crucial physics involved in the 
fragmentation of such shocked, stream-fed gas. 

These cosmic filaments play a pivotal role in galaxy 
formation by adding a third accretion mode to the two 
processes that have traditionally been thought to par- 
take in galaxy formation. Spheroidal galaxy components 
were thought to have been built by gas channelin g in 
galactic mergers ( Toomre|1977 White fc Rees|1978| and 
disk components by accretion from a shock-he ated in 



a first galaxy to be one of the first DM halos to form 
with mass > 5 x lO 7 A/0, thus meeting the requirement 
that gas heated to the virial temperature is able to cool 
through hydrogen excitation (or Lyman-a) emission. 

A further motivation for this work is understanding 
the formation of low-luminosity dwarf spheriodal satel- 
lite galaixes (dSph) and globular clusters (GCs). These 
structures have similar luminosities, stellar masses, gen- 
erally no gas nor recent star formation, and contain 
metal-poor, old stellar populations. However, GCs are 
small (~ 1 — 10 pc) and do not seem to contain much 
DM while dSph galaxies are much larger (> 100 pc) and 
are heavily dominate d by DM as indicated by their larg e 



mass-to-light ratios (Mateo 19981 Strigari et al 



It is clear these are two structurally distinct p opu 



2008 



at ions 



(e.g., Kormendyl 1985 Belokurov et al. 2007), but it is 



tergalactic medium (IGM ) during virialization (Rees & 
Ostriker|1977 Silk|1997 1 . Hot accretion occurs when gas 



is shock heated to approximately the virial temperature 
close to the virial radius of the DM halo; the gas then 
cools, falls inward, and circularizes in the galactic disk. 
Recently, however, a new mode of accretion has been 
found to dom inate in high-cr cosmological peak halos a t 
high redshift pekel fc Birnboim|2006] |Dekel et al.|2009[ ). 
This mode, coined 'cold flow accretion, occurs along cos- 
mic filaments that feed gas directly to the inner parts of 



necessary to understand the cosmological circumstances 
of their formation. To answer whether the first galaxies 
are progenitors to dSph galaxies, GCs, or both requires 
an understanding of star formation in these objects. 

With this in mind, we investigate the fragmentation 
properties of the metal-pre-enriched cold flow accretion 
in protogalaxies at high redshift, and explore implica- 
tions for star formation in these objects. Dense, metal- 
enriched baryonic streams flow along the filaments of 
the cosmic web and penetrate deep into a protogalaxy. 
The multiple (e.g., ternary) streams collide with each 
other, or collide with a turbulently-supported gas accu- 
mulation that has already collected in the center of the 
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protogalaxy. Since the stream inflow is supersonic, the 
colliding gas must pass a strong shock, but then cools 
rapidly under approximate isobaric conditions. We iden- 
tify the characteristic fragmentation mass scale in the 
shocked, isobarically cooling gas. We find that if a strong 
photodissociating Lyman- Werner (LW) background has 
significantly suppressed molecule formation, the charac- 
teristic mass scale depends strongly on the metallicity 
content of the gas. 

The outline of this work is as follows. In Section [2] 
we describe our methodology for the cooling of shocked, 
metal-enriched gas. Section [3] motivates the initial con- 
ditions for our model, including metal abundances and 
the physics of cold accretion. In Section [4] we discuss our 
fragmentation criteria and the results from our numeri- 
cal integrations, and summarize our findings and broader 
implications in Section [5j Throughout the paper we as- 
sume the ACDM cosmological parameters: fl m = 0.27, 
Oa = 0.73, an d h = 0.70, consiste nt with the 5 yr. 
WMAP results ( |Komatsu etaL|2009 1 . 

2. COOLING OF SHOCKED, METAL-ENRICHED GAS 

A cold accretion stream penetrates deep inside the pro- 
togalactic halo until it encounters an accretion shock. 
Subsequently, the sh ocked gas layer grows in size while 



cooling isoba rically (Shapiro & Kang 1987 Clarke & 



Bromm 20031, until it becomes gravitationally unstable 



and fragments into Jeans-mass-sized clumps. Our objec 
tive is to calculate the initial masses of the fragments and 
we do not follow their evolution past the point of initial 
fragmentation. To this end, we model the evolution of 
the shocked, cooling gas accurately within the one-zone 
approximation. We utilize a comprehensive chemical net- 
work to track the time evolution of the abundances of all 
primordial species (H, H + , H , H2, H^", e~, He, He + , 
D, D + , and HD) and the most important metal coolants 
expected at hi gh redshifts: neutral and singly ionized O, 
Si, C, and Fe dBromm & Loeb||2003a; Santoro & Shull 



2006[ |Maio et al.|p07[ |Glover fc Jappsen||2007[ ). TKe 

cooling gas is subject to compression by the ram pres- 
sure of the accretion shock which suggests isobaric con- 
ditions are app roxi mat ely valid in the post-shock region. 
In Sections |2 . 1 1 and [2 . 2| we proceed to discuss cooli ng a nd 
heating processes in the shocked gas. In Section [273] we 
discuss the relevant ch emic al processes in the post-shock 
gas. Finally, in Section 2.4 we describe the time evolution 
of the shock-compressed gas under isobaric conditions. 

2.1. Cooling Processes 

We calculate all relevant cooling and heating processes 
in low-density, shock heated gas at high redshift. The 
cooling of gas with primordial abundances set by big 
bang nucleosynthesis has been extensivel y studied (e.g., 
Shapiro fc Kang|l987| [Oh fc Haiman|2002| ) . Above 10 4 K, 
hydrogen and helium excitation are the dominant cooling 
channels. Between 200 K < T < 10 4 K, if molecular hy- 
drogen forms, the cooling that it facilitates is dominant. 
If present, HD can continue to cool the gas below 200 K 
due to its permanent dipole moment and smaller rota- 
tional energy level spacing, although an elevated electron 



fraction is necessary for its efficient formation flUehara 
fc Inut suka 2000; Naka nuTra fc Umem ural 
fc Bromm||2006[ |Yoshida et al.||2007 ) 



oimson 
etals, supplied 

trom previous star formation, can be a significant coolant 



at T < 10 4 K due to forbidden, fine structure photon 
emission from low-lying electronic radiative transi tions 



( [Sutherland fc Dopita|19"93}|Gnat fc Sternb erg 200 7} and 
can substantially change the gas dynamics ( |Omukai|2000| 
Bromm & Loeb 2003a). Below 10 4 K, our cooling rate 
as a function of metallicity closely matches the cooling 
rates calculated by other au thors in gas of similar tem- 
peratures and densities (e.g., Santoro fc Shull|2006 Maio 



et al. 



20071 



Smith et al. 



2008) 



ture emission from O, Si, Fe 



We consider tine struc 
and C — the species that 



dominate metal cooling in a low density, low metallicity 
gas (e.g., ISantoro fc Shull||2006| ). We find that in the 
shocked gas m an atomic cooling halo, Si, Fe . an d C are, 
to a large extent, singly ionized (see Section 3.3), but O 



remains mostly neutral. Above 10 4 K, our cooling rate is 
not accurate at high metallicities (especially for Z > Zq) 
as we do not include the multiply ionized species, but 
since we restrict ourselves to Z < 10~ 2 Z (see Section 
k3J) , t his assumption does not affect our results (see Figure 



Gnat fc Sternberg 2007) 



Due to IN LI E conditions when simulating galactic ac- 
cretion from the IGM, Boltzmann statistics does not de- 
scribe the state of an atom or molecule, but instead, 
level populations have to be calculated from the prin - 
ciple of detailed balance (see Glover & Jappsen 2007), 
where electronic collisional excitations and de-excitations 
are balanced by spontaneous downward transitions. Ex- 
citation and de-excitation rates, energy level spacings, 
and spontaneous emission coefficients for C, C + Si, Si + , 
Fe, Fe + , O, and 0+ are the parameters needed to deter- 
mine the fine structure co oling rate. We refer the reader 
to other sources, such as [H ollcnbach fc McKeel (|198 9[ 
(hereafte r HM8 9) , |Santoro fc Shull| d2U0(j| ) , and |Glover 
Jappsen] (2007), for detailed discussion of the individual 
cooling rates. We here only provide the refere nces from 
which we co mpiled our data: C - HM89 and |Johnson| 



et al.| (|l987|; C 4 



si^ 



Glover & Jappsen (2007) and" Bell et al.| 
HM89; IRoueffl (|1990|); Fe, Fe + 



HM89 and Wilson & B elli fl2002| ); O 
|TggS| ; Si 
- HM89. 



HM89; si 

The simplest model of fine structure cooling from an 
atom is one in which the electron has only one excited 
state to occupy, the so-called two level system. In this 
syste m, the excited st ate level population is given by 
(e.g., |Maio et al.||2007[ ) 



"2 



7l2 +7l2"c/^H 



7? 2 



721 + (7i2 + 72iH/«H + A 21 /n H 



"tot , (1) 



where n 2 is the number density of particles in the sec- 
ond excited state, n to t the total number density of the 

species, while 7i2* an d 72i I ' e ^ are > respectively, the ex- 
citation and de-excitation rates due to collisions with 
neutral hydrogen and free electrons (collisions with hy- 
drogen ions, H2, and other species are subdominant). 
The corresponding volumetric cooling rate for fine struc- 
ture line emission is given by Aii no = r^^iAE^i, where 
A21 is the spontaneous emission coefficient and AE 2 i is 
the energy difference between states. For n level sys- 
tems, n coupled equations are solved to determine the 
level populations. The CMB temperature sets a t hermo- 
dyna mic lower limit on the gas temperature (see Larsonj 
19981; we incorporate this limit by replacing the cooling 
rates with the effective rate: A off (T) = A(T) - A(Xcmb)- 
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Fig. 1. — Post-shock non-equilibrium cooling rate for differing 
metallicites for a 10 8 Mq at z = 10, using a solar pattern of met- 
als. This figure only includes cooling from hydrogen excitation and 
metallic fine structure emission - molecules are ignored. The bump 
at T ~ 200 K is due to the density increase towards lower temper- 
atures due to our isobaric density evolution. The initial density 
is n ~ 4 X 10 3 cm -3 . Towards 30 K the cooling rate drops off as 
we approach Tqmb- Below T ~ 10 4 K the cooling rate is a strong 
function of metallicity. 




100 1000 10000 

T [K] 

Fig. 2. — Post-shock non-equilibrium cooling rate showing the 
contributing components for an accretion shock in a 10 s Mq DM 
halo at 2 = 10 with J21 = 10 and metallicity Z = 10 — 2 Zq with a 
solar pattern of metals. It is evident that fine structure emission 
from Fe + , Si + and O are the dominant metal coolants. C + is 
not effective as the density is above its critical density, n cr c + ~ 

3xl0 3 cm — 3 . H2Cooling dominates from 1000— 6000 K while metals 
dominate below 1000 K. HD is always an insignificant coolant at 

io- 2 z Q . 



For illustration, in Figure [T] we show the non- 
equilibrium cooling rate as a function of metallicity, eval- 
uated at every step during the evolution of the post-shock 
gas in our fiducial model with a 10 8 Mq halo at z = 10, 
solar distributions of metals, and with a h igh L W radia- 
tion background ( J21,lw = 10 4 , see Section 3.3 1. The lat- 
ter significantly photodissociates H2, and will keep metal- 
lic species singly ionized (see Section 3.3 ). The wide peak 
at T ~ 200 K is due to the isobaric increase in density 
towards lower temperatures. The initial values for den- 
sity and temperature are 4 x 10 3 cm -3 and 1.1 x 10 4 K, 
respectively, and the initial chemical abundances are set 
to the values corresponding to thermal and statistical 
equilibrium at this temperature. Below T ~ 10 4 K, the 
cooling rate is a strong function of metallicity. Fig- 
ure[2]shows the individual elemental contributions to the 
nonequilibrium cooling curve for a solar abundance pat- 
tern. Metal cooling, particularly from O, Fe + , and Si + , 
dominates below 800 K; H2 dominates between 800 and 
8,000K .The critical density for the C + fine structure 



transition is n 



cr.CT 



A10/7: 



in 



3 x 10 3 cm" 



Thus, 



at the densities sampled by our isobarically cooling gas, 
n = 10 3 — 10 6 cm -3 , the C + line is in LTE, which makes 
it a weaker coolant compared to Fe + , Si + , and O, whose 



critical densities are 



10° 



10 cm 3 so thermaliza- 



tion o ccurs at a later evolutionary stage ( Santoro & Shull 
2006 1. In this regime, HD is always an ineffective coolant 



If, however, metallicity is Z < 1O _5 Z0, and molecules 
are permitted to form, the gas is able to reach Tcmb by 
HD cooling. 

2.2. Heating Processes 

Since we here ignore the effect of dust, the predomi- 
nant source of heating is UV pumping of H 2 , particularly 
at higher densities. To compute the rate, we assume a 
pumping rate 8.5 times that of the photodissociation of 



H 2 ( |Draine fc Bertoldi|[l996| |Glover fe Jappsen||2007 ). 
If each pumping pr oduces ~ 2 (1 + n CT /n)~ L eV of heat 
( Burton et al.|1990 ), where n cr is the critical density for 
H2 (which we always take to be 6, 000 cm -3 ), the corre- 
sponding heating rate is: T = 2.7 x 10 -11 i?diss^H 2 (1 + 
r/n)^ 1 ergs -1 c m~ 3 ; the computation of i?diss is dis- 



cussed in Section 13731 below. 

Other significant sources of heating are the photoion- 
ization of metallic species and the radiative coupling to 
the CMB. It is safe to assume that at z ~ 10, prior to 
complete reionization, the mean free path of H ionizing 
photons is very short compared to photons below the Ly- 
man limit, thus rendering heating from H and He ioniza- 
tion unimportant. However, photons below the Lyman 
limit should be l // eS c ~ 10 2 times more abundant than 
H-ionizing ones ( |Bromm fc Loeb||2003b |, making pho- 
toionization heating of metallic species with first ioniza- 
tion potentials below 13.6 eV (C, Si, and Fe) dominant. 
This also has consequences for the photodissociation of 
H2, whi ch o ccurs in the LW bands just belo w 13.6 eV (see 



Section 3.3). As discussed in Section 2.1 CMB heating 
is achieved by modifying the effective cooling rate. Over- 
all, heating of the gas is generally not a significant effect 
for the evolution considered here. 

2.3. Chemistry 

Our chemical model includes all relevant primordial 
chemical species: H, H+, H- H 2 , H+ e" He, He+, 
Hc++ 



D, D + , and HD (Haiman et al 



Palla| |2002 1 

coolants- 



1996 



Galli & 



We also include the most important metal 



neutral and singly ionized C, Si, O, and Fe. 
Collisional ionization, recombination, and charge ex- 
change reactio ns for C, Si, and O are taken from the 
compilation in Jappsen et al. (|2009b[). The reactions 



for Fe were taken fro m Hollenbach fc McKee (1989) and 
Woods et al. (1981). We reproduce the complete cat- 
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alog of rates in our chemical network in A ppendix A 
usi ng the s ame numbering convention as in Jappscn et 



al 



(|2009bj for rates 1-222; after 222, we start our own 
numbering convention. We initialize our chemical model 
with equilibrium abundances atT = l.lxl0 4 K. These 
initial abundances are njj - 4x 10 3 cm~ 3 , n H + /njj — 35, 
n H -/nn — 3 x 10 -5 , n H2 /n H ^ 10~ 2 , n H +/n H ^ 10~ 5 , 

n e /n H ^ n H +/n H ^ 30, n H D/«-H ^ 10~ 6 , n Hc +/n H ^ 
8 x 10~ 3 , and n He ++ /rin - 4x 10 -9 . Metal species are 
initialized such that Si, Fe, and C are fully ionized while 
O has the same ionization fraction as hydrogen. 

Chemically, the formation and destruction of H 2 /HD 
has the strongest impact on our results as it can be a 
very effective low-temperature coolant. In the absence of 
dust, molecular hydrogen can form through t wo channels 
(e.g., |Abel et al.|[l997[|Tegmark et al.||l997[ ). In the H" 
channel, H~ and H combine to form H2 and e~ , whereas 
in the channel, and H participate in a charge ex- 
change reaction producing H2 and H + . The efficiency of 
these formation channels strongly depends on the abun- 
dance of f ree electrons that are required for H _ and 
synthesis (Oh & Haiman 2002). In all our calculations, 
the H channel dominates. Destruction of H2 is mainly 
from UV phot ons in the LW bands and will be discussed 
in Sec tion [3J3I HP may be an important coolant below 
200K (Johnson & Bromm 2006). Its formation proceeds 
by an exchange reaction between D + and H2. In our 
physical scenario, HD will only be a significant coolant 
in a low metallicity environment. 

2.4. Isobaric Cooling and Time Evolution 

If the post-shock cooling is thermally stable and if the 
shocked gas is in approximate pressure balance, the ther- 
modynamic evolution of a gas packet that has passed the 
shoc k obeys conservation o f energy, as expressed by (see, 
e.g., |Shapiro fe Kang||1987[ ): 



d(pe) 
dt 



p\ dp 
dt 



A + F, 



(2) 



where p, e, p, A, and T, are, respectively, the gas den- 
sity, specific internal energy, pressure, volumetric cool- 
ing and heating rate. To express these in terms of the 
chemical species number densities ni, adiabatic indices 
7i, and the temperature T, we write: p = ^ nik^T and 
pe = J^. Uik^T I (7j — 1). The density derivative dhip/dt 
can be written as the total number density derivative 
dlnn/dt, where n — X)i n ii minus the number density 
derivative due to pure chemistry evolution (isochoric), 
as follows 



d In p 



d In n 



dt \ dt 
d In n 



dt 

d In T 

dt 



fixed 

d In n 
dt 

d In n 

dt 



ochori( 



isochoric 



(3) 



where in the last step, we imposed isobaricity, 
din n/dt = —din T/dt. 'Fixed' refers to evolution with 
fixed abundances. With this, equation ^ becomes 



E 



linjk-Q 
7,-1 



dt 



- A + r (4) 

We further calculate the change of species abundances 
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V ^ / isochoric- 
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dt 



1, 



,iV s1 



.)■ 



isochoric 



(5) 

Equation Q and (pi) are solved for dT/dt and dnt/dt, 
and the resulting ^v^p CC ios + 1 differential equations are 
integrated implicitly using the Bulirsch-Stoer-t ype, semi 



Bader & 



lme step 



implicit extrapo lation mid-point method (see 
|Deuflhard||1983[ ) ■ The coarsest Bulirsch-Stoer 
is set to At = min{ni/|ni|}. 

Isobaricity, i.e., pressure equilibrium, is a valid approx- 
imation if the sound crossing time of the post shock 
region is shorter than the cooling time, £ SOU nd < ^cooi- 
In Figure [3j we confirm that this condition is satisfied 
throughout the evolution of the gas leading to grav- 
itational instability and fragmentation, except for the 
highest metallicity case, Z = 10~ 2 Z Q , when the cool- 
ing time becomes shorter than the sound crossing time 
at ~ 2 x 10 4 cm~ 3 . We thus conclude that for high 
metallicity, thermal instabilities may play a dynamically 
significant role. It is interesting that these instabili- 
ties only occur in the prese nce of efficient metal co oling. 
Thermal instabilities (e.g., )Fieldfl965| |Balbus||1986| are 
one proposed mechanism for the formati on ot molecu- 
lar clo uds and the origin of turbulence. |Niklaus et al.| 
(20091 performed two-dimensional adapative mesh re 



finement (AMR) simulations of colliding flows with non- 
zero metallicity and found turbulent, clumpy structures 
formed from thermal instabilities triggered by a super- 
sonic collision. It is tempting to suggest that the thermal 
instabilities we anticipate here at substantial metallicities 
could lead to dumpi ness and turbulence like that seen in 
Niklaus et all d2009b. 



3. MODEL PARAMETERS 

The choice of parameters for our modeling of the accre- 
tion flow into the protogalaxy will govern the resulting 
fragmentation mass scale. This section focuses on the 



three most important model ingredients. In Section 3.1 
we discuss the properties of th e cold flows into high red- 
shift galaxies. In Section |3.2[ we describe the possible 
abundance patterns of accreting gas and place an upper 
limit on the metallicity content of the gas. Finally, in 
Section [3~3} we discuss the effect of a LW radiation back- 
ground on the dissociation of molecular hydrogen and 
the ionization of metallic species. 



3.1. 



Birnboim & Dekel 
virialized gas from the 1G 



Cold Flow Accretion 
2003 



consider accretion of non- 
into a DM halo and find 



that, at redshifts z > 2, under the assumptions of spher- 
ical symmetry, a stable virial accretion shock does not 
exist in halos smaller than 10 11 M due to the rapid 
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Fig. 3. — Test of the validity of our isobaric density evolution 
assumption. The four lines correspond as shown to the four cases 
of Table [2] As shown, isobaricity holds except in case A (strong 
metal cooling), where it is violated for densities above ~ 10 4 cm -3 . 
This will likely lead to runaway cooling collapse in the density 
regime where t CO olAsound < 1. In the other cases, where H2 cooling 
dominates, we conclude that pressure support acts on a smaller 
timescale than radiative cooling. This effectively allows the ram 
pressure of the cold flow to keep the shocked region at approxi- 
mately constant pressure throughout the evolution, at least until 
fragmentation occurs. Runs C and D undergo fragmentation be- 
fore being able to reach high densities, and we do not follow their 
evolution past this point. 



cooling of the post-shock gas. Our model halos are well 
within this regime. Instead, the accreting gas contin- 
ues on supersonic infall traject ories until it reaches much 
closer to the center of the halo 
argue that such 



Dekel k Birnboim| ([2006 ) 
'cold streams provide a natural expla- 
nation for many observed galactic phenomena, includ- 
ing the existence of starbursts lacking the evidence of 
a major merger, the existence of bulgeless disks, and 
the blue-red bimodality observed in present-day galaxies. 
Keres et al. (2005) find, using smoothed particle hydro- 
dynamics (SFH) simulations, that in halos with masses 
< 10 Mq at z = 3, almost none of the gas that enters 
the galaxy has been heated to the virial temperature, 
further supporting the conclusion that cold accretion is 
the main feeder of baryonic material into galaxies at high 
redshift. 



Wise k Abel ( 2007 ) simulated the cosmological forma- 
tion of a ~ 10 8 M & cosmological object and find that 
the virialization of shocked accretion streams creates a 
nearly isothermal (~ r~ 12 / 5 ) density profile co ntaining 



Greif et al. 



(|2008| carried 
10 and found 



supersonically turbulent gas 
out a simulation a similar mass halo at z 
that filamentary accretion dominates the buildup of the 
baryonic component of the galaxy. In their simulation, 
the density of cold accretion streams just prior to im- 
pact with the central turbulent mass concentration was 
~ 10 3 cm~ 3 with a temperature of ~ 200 K. The pre- 
shock number density in cold streams can be estimated 
by considering infall in a halo with mass M — 10 8 Ms Mq 

and virial radius i? v i r = 1.5Mg^ 3 [(l + z)/10] _1 kpc, 
where the latter we define as the radius within which 
the mean density of a halo equals 200 times the critical 
density of the universe. If the density of the infalling 



extrapolated to the innermost 10 pc, the approxima te ra- 
dius of cold-stream collisions and shock formation ( Wise 



et al |2008| ), is 
result. 



10 cm , consistent with the numerical 



3.2. Abundances from Early Enrichment 
We restrict ourselves to an upper metallicity limit of 



10 2 Zp, . This follows from the work of Greif et al. 



(2008) 



who find that a 10 8 Mq halo will be assembled out of ma- 
terial that^lias_e2crjeriejiced a maximum of ~ 10 Pop III 
SNe. |Bromm et at] |2003[ ) showed that in a Pop III SN 
explosion, ~ 907o~of the metall ic yield will be ejected 
into the IGM. Greif et al. (2010) find that accretion in- 
flows return roughly 50% of the metallic enrichment to 
the central halo by the time it grows to 10 s Mq . Assum- 
ing a homogeneous density distribution of metals over 
the halo, roughly 1 proper kpc in size, gives a metallicity 
fraction of ~ 10~ 3 Zq, suggesting that 10 -2 Zq is a con- 
servative upper limit for the metallicity. We explore the 
effect of metallicity that varies from a primordial compo- 
sition, which is effectively described by Z < 10 -8 Zq, to 
this upper limit, using both PISN and solar abundance 
patterns. 

There have been assertions that a LW background from 
Pop III stars could inhibit further star formation by de- 
stroying H 2 and H~, which would keep the first galax- 
ies starless and thus metal free or metal d eficient (e.g., 
Haiman k Bryan|2006| |Dijkstra et al.|2008[ ). Other stud- 
ies have shown that an ionizing background enhances 
metal enrichment as a hig h rate of ionization en courages 
molecule formation (e.g. 



Johnson et al. 20081 and en 



hanced cooling. Another recent study revisited the ques- 
tion of radiative feedback from Pop III stars using de- 
tailed radiation-hydrodynamical simulations and found 
radiative feedback to be neutral (Ahn & Shapiro 20071 
for Pop III. 2 star formation (stars still forming from pure 
H/He gas but with radiative feedback from previous star 
formation), suggesting there will be a metal abundance 
of at least 10 Zq in these protogalaxies, depending on 
the halo capture fraction of SN nucleosynthetic prod- 
ucts. The metal a bundance pattern i s also not k nown 
with any certainty. Umeda k Nomoto ( 2002 2003 1 sug- 
gest that the abundance pattern reflects enrichment from 
hypernovae, produced by stars in the 20 — 120 Mq mass 
range, while other authors argue that Pop III stars end 
their lives as PISNe, producing a very different metal- 
licity pattern (Heger & Woosley 2002). Therefore, it is 
crucial to explore a range of abundance patterns in our 
modeling. 

Furthermore, the mixing of heavy elements within the 
protogalaxy after they get incorporated through cold ac- 
cretion streams, or following a SN explosion inside the 
galaxy itself, remains poorly understood. Heavy ele- 
ments produced from Pop III stars may never return to 
the galaxy and may reside in the IGM indefinitely. An- 
other possibility is that metals may leave the host galaxy 
after a SN and return at a later time. Once the galaxy 
is assembled, will there be a region of high metallicity 
contained around SNe, or will metallicity quickly dif- 
fuse throughout the protogalaxy? Inhomogeneous mix- 
ing may create pockets of enhanced enrichment, with a 
metallicity much greater than Zq . The degree of mixing 
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or clumpiness of the metal distribution inside the proto- 
galaxy will be investigated in greater detail in later work. 

Our nu merical integrations are run for both standard 
solar (see Grevesse & Sauval 1998) and PI SN metal 
abundan c e patt erns. PISiN patterns are from |Heger fc] 
Woosley| p002|) for a 150 and 250 M Q progenitor SN 



s corresponding to a 75 and 125 M helium core, respec- 
tively). Table Vu represents the number density (in cm~ 3 ) 
of metallic species given a hydrogen number density of 
lcm -3 . The PISN values are computed with the same 
overall mass fraction in metals, but differing elemental 
proportions compared to the Sun. 

3.3. Radiation Field 

The rate of photoionization of a s pecies (-Bphoto) fo r an 
optically thin gas is given by (e.g.,|Osterbrock||1989[) 



R 



photo 



hv 



o~(v) dv, 



(6) 



where J v = 10 -21 J 2 i ergs s _1 cm -2 Hz~ sr _1 is the 
mean intensity, o(v) is the frequency dependent pho- 
toionization cross section, and the integration runs from 
the threshold photoionization frequency, t'th, to infin- 
ity. In the case of me tallic photoionizations which are 
dominant (see Section 2.2 1, we run this integration to 
vh, the ionization threshold of hydrogen, as hydrogen 
ionizing photons are effectively shielded by the IGM 
prior to reionization. Photoionization cross-s ections for 
C, Si, Fe , and O are taken from the fits of |Verner et 
al. (19961. For our incident spectrum, we use thermal 
Planck spectra with T* = 10 4 and 10 5 K and a powe r 
law spectrum, J v oc v~ x . As in Bromm & Loeb (2003a 
we normalize these spectra at the Lyman limit by J 2 
The photodissociation of H 2 is a more complicated pro- 
cess. It occurs through t he two-step Solomon process 
( Stecher fc Williams||1967 ) by photons in the LW bands, 
'. 1.2— 13.6 eV. A commonly used approximation for this 
rate of dissociation, -Rdiss,H 2 i is given by ~ 1.38 x 10 9 J P 
where J p is the mean intensity a t 12.4 eV, the mea n en- 
ergy of a LW band photon (see Abel et al. 1997). We 
assume the photodissociation rate ot HD equals that of 
H 2 . 

At high column densities, H2 can be particularly ef- 
fective at shielding itself from LW radiation whi ch can 
significantly a lt er its photodissocia tion rate (e.g., Oh & 
|Haiman|[2002 l. |Draine fc Bertoldil ( p96| ) provide an ap- 
proximation for the shielding factor such that i?diss oc 
Ahwd, where / shicld ~ min[1.0, (7V H2 /10 14 c m - 2 )- - 75 ] 
with 7Vh 2 being the column density of H2. However, 
this applies to a static, cold medium. If the LW bands 
are Doppler shifted by large thermal or bulk motions, 
the shielding rate will be much lower. In in accreting 
protogalaxy, where large-scale velocity gradients reach 
~ 20 km s _1 , self-shielding will certainly be overesti- 
mated within this pres cription. Taking thermal b roaden- 
ing into consideration, Draine & Bertoldi ( 1996 ) provide 
another fit to the shielding factor as 



,/shield — 



0.965 



0.035 



7- ' 77 " exp[-8.5xl0- 4 (l+x) - 5 ] 

(l + x/b 5 ) 2 (l + x)°- 5 ^ L 



where x = A^ H , and b 5 = b/10 5 cm s 1 . We follow Ahn & 



ill 



Shapiro (2007) in setting b = 9.12 km s 1 , as their phys 



ical system is approximately similar to ours. Although 
this will still overestimate the amount of self-shielding 
when bulk motions dominate, it is an improvement over 
the previous fit and will be used in computing our value 
of H 2 self-shielding. 

The radiation field present in these high rcdshift pro- 
togalaxies is uncertain, owing to the contribution from 
the external, global IGM field, and the local field pro- 
duced by internal Pop III/II stars. One can estimate, 
at z = 10 , that the intensity just above the Lyman 
limit, Jjj, necessary to reionize the universe is ~ 40 
(Bromm & Loeb 2003a). However, photons with ener- 
gies below the Lyman limit are virtually free to escape 
their host halo, so that the corresponding intensity, J^, 
which includes the H 2 dissociating LW bands, could po- 
tentially be J" 1 ~ 100 times higher, w here / e sc is the 



escap e fraction of H-ionizing p hotons (|Wood fc Loeb 
20001). Other recent estimate s ( |Razoumov fc Sommer 
Larsen|2009||Wise & Cen|2009|) have found f RSC , ~ 0.25- \ 



in dwarf galaxies at high redshift. Finally, the LW back 



global aven 
|aH [20081 IW 


ige LW field in clustered halos ( 


Dijkstra et 


halen et al. 


2008; Johnson et al. 


2008; |Ahn 


et aT|2009 


Shang et al. 


2009 


) . For this study, we take 



J21 = 10 and J 2 i = 10 4 , representing small and large ra- 
diation fields which could be realistically present at this 
epoch. 

4. THERMAL EVOLUTION AND FRAGMENTATION 

We now consider the case of a cold accretion stream 
which is shock heated in t he center of a primordial pro- 
togalaxy. In Section |4.1| we derive a simple criterion 
for when the post-shock region will fragment based on 
the free-fall and sound crossing times. Section |4.2| de- 
scribes the thermal evolution of the post-shock region, 
and presents the fragmentation m ass s cale for different 
parameter choices, whereas Section [4. 3| outlines a simple 
analytical model for the critical metaliicity in molecule- 
free regime. 

4.1. Fragmentation Criterion and Mass Scale 

To derive a rough estimate of when fragmentation will 
occur, we consider the growth of the shocked gas region 
within our one-zone model. Following the shock that 
terminates the accretion stream near the center of the 
protogalaxy, the gas begins to cool off isobari cally (see 
Shapiro fc Kang||1987| |Clarke fc Bromm||2003[ ) from the 



virial temperature of the halo, T v ; 



2 x W 4 n , 6 M% /3 [(l + z)/10] K, where \i = 0.6 /i .e is the 
mean molecular weight in the immediate post-shock gas. 
Applying the standard Jeans analysis for gravitational 
instability, we take it that fragmentation occurs when the 
ever-increasing sound crossing time equals the decreas- 
ing free fall time of the post-shock region. Specifically, 
we examine three time scales: the cooling time, t cao \ ~ 
§/c B Tn t ot/A(n,T), the free-fall time, tg = (Sn/WGp) 1 / 2 , 
and the sound-crossing time, t s ~ L/c s , where L is the 
length of the shocked slab, A(n, T) the volumetric cool- 
ing rate, and c s = {^fk^T '/ '/inij,) 1 / 2 the sound speed. The 
pre-shock density po is related to the instantaneous post- 
shock density p 1 via po(v in{aa + v shock ) = piw sh ock- Since 
the shock is strong, pi ~ 4p - If the post shock gas 
further cools and isobarically compresses to p(t), then 
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TABLE 1 

Metallicity Yield Patterns. These values represent, for a gas with a 

hydrogen number density of lcm" 3 , the number density of metallic 
species wtih differing abundance patterns normalized to the metallic 
mass fraction of the sun. 



Pattern 


C 




Fc 


O 




Si 




Solar 


2.5 x 10- 


-4 


2.8 x 10~ 5 


4.6 x 10 


-4 


3.2 x 10" 


5 


150 M PISN 


5.0 x 10" 


-5 


1.2 x 10" 6 


5.1 x 10 


-4 


1.4 x 10" 


-4 


250 M PISN 


1.4 x 10" 


-5 


1.3 x 10~ 4 


1.5 x 10 


-4 


1.0 x 10- 


-4 



TABLE 2 

Parameters for the Different Runs 



Run 


Metallicity 


J21 


T* 


M B e Final 


A 


10- 2 Zq 


10 4 


10 4 


~3M 


B 


10- 6 Zq 


10 


10 4 


~ 10 M 


C 


1O" 6 Z 


10 3 


10 4 


~ 10 5 Mg 


D 


10- 6 Zq 


10 3 


10 s 


~ 10 4 M Q 



the shock velocity in the frame of the gas will be fur- 
ther reduced by a factor of ~ p\/p{t), or u s hock/finfaii ~ 
pap\/\p{t)(px - pa)} ~ jp /p(t), and thus the length of 
the post shock region is (ignoring numerical factors order 
unity), 

L(t) ~ v shock t ~ -~ v vir t, (8) 



23M R 1/3 [(1 



where w in faii = v y - a = (GM/R^) 1 / 2 
z)/10] 1 / 2 km s _1 the virial (or infall) velocity. Isobaric 
cooling requires i coo i > t s , while gravitational (Jeans) in- 
stability and the onset of fragmentation requires t s > 
It is possible that gravitational instability does not set 
in before the gas cools to the CMB temperature. In 
this case, isobaricity is guaranteed and the fragmenta- 
tion mass scale is entirely set by the CMB temperature 
and the initial conditions of the accretion flow. 

We assume that following the initial gravitational in- 
stability, the gas fragments into clumps with masses 
given by the local Bonnor-Ebert (BE) mass, 

A/be ^ 150 M Q (T/200K) 2 , (9) 

assuming isobaric density evolution with Ti = T vrr = 
1.1 x 10 4 K and n t = 4.0 x 10 3 cm . The resulting clump 
size sets an upper limit to the masses of the star in the 
aftermath of the fragmentation. Due to subtle processes 
affecting post-fragmentation accretion, the final stellar 
masses could be a fraction of the fragmentation scale: 
M* ~ oMbe, where a < 0.5 in s tar formation in the 
local universe flMcKee k Tan||2002[ ) . 



4.2. Case Studies 

For our case studies, we consider a high and low metal- 
licity case (10 -2 and 10 -6 Zq), a high and low LW ra- 
diation fields (J21 = 10, 10 3 , and 10 4 ), and source tem- 
peratures of T* = 10 4 and 10 5 K. We do not present a 
run with high metallicity and efficient molecular cooling, 
but we do comment on the outcome of such a run at 
the end of this section. In Table [2] we provide the pa- 
rameters of our four representative runs. We present the 
results for each case in succession. See Figures [4^[7] for 
the temperature history, cooling rate, H2 fraction, and 



ionization degree (defined such that ionization degree = 
n e - free /n p + any form), respectively, for each run.. 

Run A represents a relatively high metallicity environ- 
ment with a strong LW background that photodissoci- 
ates H2, rendering metals the predominant coolant at 
T < 10 4 K. As in all our runs, we assume that the 
baryons in DM filaments penetrate deeply into the center 
of a protogalaxy where they undergo a strong shoc k with 
M. = fvir/cs ~ 10. Calibrating to th e gas density in Greif 



et al. (2008) as explained in Section 3.1 above, the post 
shock number density is 4no ~4x 10 3 cm~ 3 . The shock 
raises the gas temperature to T v ; r ~ 1.1 x 10 4 K. At this 
temperature, the gas cools effectively via atomic hydro- 
gen line excitation, while the size of the post-shock region 
begins to increase in size according to equation After 
~ 10 5 yr the gas has cooled to ~ 35 K, which is slightly 
above Tqmb(z = 10) ~ 30 K due to UV pumping of H 2 , 
and after ~ 1.5 x 10 6 yr the post-shock region has reached 
the fragmentation criterion (tg < t SOU nd)- We do not 
follow the evolution past this point. The radiation back- 
ground (J21 = 10 4 ) has prevented an appreciable amount 
of H2 from building up, with /h 2 = riH 2 / ri H never exceed- 
ing 10~ 7 . The BE mass at the point of fragmentation is 
~ 3M Q . Using a fit to the baryonic growth r ate of ha- 
6.6M& 5 (1 + z) 2 - 25 /.i 65 M yr- 1 jDekel et 
20091), where M 12 



los of M 



M vir /1O 12 M and f. 165 is the 



o baryonic fraction in units of the cosmological value 
^b/^m — 0.165, we find that the protogalaxy accreted 
~ 5 x 10 4 Mq of baryons before the onset of fragmenta- 
tion. If all these baryons experienced fragmentation, this 
would represent a region with ~ 10 4 separate fragments, 



Bromm k Clarke 


20021. For more 


tion efficiencies, 


t order ~ 0.1, wc 



sizable cluster of Pop II stars to arise. 

Run B explores a low metallicity environment together 
with a feeble LW background which allows for H2/HD 
formation and cooling. Starting from the same initial 
density and temperature as Run A, molecular hydrogen, 
now able to form more effectively, reaches an abundance 
/h 2 ~ 10~ 3 only 10 3 yr after the shock passage; H 2 and 
HD dominate the cooling from 8, 000 to 30 K (see Fig. [5]). 
Like in Run A, the post-shock region fragments after 
~ 10 6 yr at a temperature of ~ 50 K. The final BE mass 
is ~ 10 Mq. Overall, this run is very similar to case A, 
but cooling is now dominated by H2 and HD instead of 
metallic fine structure lines. The likely end product is a 
collection ~ 10 3-4 fragments. 

Run C is intended to illustrate the result for a post- 
shock region with little metal enrichment and a large 
enough LW background to kee p H2 cooling suppressed , 
similar to the case examined in Bromm & Loeb (2003b). 
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Fig. 4. — Post-shock temperature history of four different runs, 
all performed with a solar abundance pattern. See Table[2]for a de- 
scription of the run parameters. All runs start at T = 1.1 X 10 4 K ~ 
T v ; r . The addition of molecules (^predominantly) produces ad- 
ditional cooling at high temperatures (T ~ 10 3 K) and thus the 
temperature falls off from T ~ 6500 K first. Metals, which are ef- 
fective low temperature coolants, allow the gas to reach the Tomb 
quickly, as shown in case A. Case C initially has no effective coolant 
below T ~ 10 4 K, but as H2 self-shielding builds up the gas is able 
to cool a small amount before fragmentation. Case D, with a black- 
body source temperature of T = 10 5 K, is able to build up a small 
Infraction and cool to ~ 500 K before becoming gravitationally 
unstable. Fragmentation occurs for each run at ~ 10 6 yr. 



Initially, the gas is still able to cool via atomic hydro- 
gen excitation to ~ 6, 000 K, when the cooling time be- 
comes similar to the Hubble time. The gas thus stalls at 
6, 000 K, at a density of 6 x 10 3 cm~ 3 , until it fragments 
10 6 yr after the shock. The final BE mass is 10 4 M while 
the total accreted mass is also of this order. This accre- 
tion episode, with no effective low temperature coolant, 
produced one large fragment, possibly indicative of an 
object which will not fragment further an d become a su- 
per massive black hole (SMBH) precurso r (|Omukai|2001| 
Bromm Ez Locb 2003a: Shan g et al.|2009[|Schleicher et al.| 
2010 \. Beca use o f our overestimation of H2 self-shielding 
(see Section |3.3| , this value reflects a lower limit on the 
fragmentation mass scale, which is most likely higher, 
possibly large enough to facilitate one large 10 5 M Q frag- 
ment capable of direct collapse to a SMBH. 

Run D has the same initial conditions as Run C, 
except that the blackbody source temperature is now 
= 10 5 K instead of 10 4 K, more represent ative of the 
type of so urces associated with Pop III stars (Bromm et 



al. 2001a). Because of the normalization at 13.6 eV, a 
blackbody of = 10 5 K is much less effective at destroy- 
ing H~ , so that H 2 formation is inhibited much less than 
for a T = 10 4 K blackbody (see FigM. The slightly more 
effective H2 cooling leads to a finalfragmentation mass 
scale of 10 3 M®. Like in case C, overestimation of H2 self 
shielding sets this as a lower limit on the fragmentation 
mass scale. 

Our results indicate that if there is a low tempera- 
ture coolant available, either metals or molecules, the gas 
will fragment with a characteristic mass of 1 — 10 M Q 
(cases A and B); a fragment in this mass range has a 
very good chance of producing a solar mass star. A run 




1000 



10000 



T [K] 



Fig. 5. — Total cooling rate as a function of temperature for the 
four representative runs. For Cases C and D, which fragment before 
reaching close to TcmBi only the pre-fragmentation cooling curve 
is shown. Metals dominate the cooling in case A at T ~ 200 K 
while H2is a very effective coolant at T ~ 3000 K. 




Time [yrs] 

Fig. 6. — Time evolution of H2, /h 2 i m the four representative 
runs. The strong photobackground in case A prevents any appre- 
ciable H2 fraction from forming, while a LW background of J21 = 10 
allows the Infraction, n(H2)/n(H), to exceed 10 — i . 



with both efficient molecular cooling and high metallic- 
ity (J21 = 10 and Z = 10~ 2 Zq) will behave similarly 
to case A and produce a ~ 1 M Q fragment. If, how- 
ever, no low-temperature coolant is available as a result 
of a high LW background and little chemical feedback 
from early SNe, there will likely be only one large frag- 
ment (~ 1O 4_5 M ). The critical LW background where 
H2 is no longer an effective coolant is approximately 
J21 = IO 3 - 4 for a blackbo dy spectrum of T * = 10 4 K, 
close to the value found by |Bromm 
similar initial conditions. | Shang et al. 
smaller critical J21, most 



Loebl (|2003b| ) for 
jQ09p find a much 
kely due to the higher H2 self- 
shielding experienced in these relatively high density ac- 
cretion flows, although as discussed, we d o consistently 
ove restim ate the degree of self-shielding. |Schleicher et| 



al. (2010) find cooling suppression for J21 ~ 1 — 10 4 , de- 



pending on the gas density. Another point of interest is 
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Fig. 7. — Ionization degree of the gas in the four runs. We define 
the ionization degree to be the number density of free electrons 
divided by the number density of protons in any form, n e - free /n p + . 



the appearance of a critical mctallicity in the absence of 
molecular coolants. In case A, the final fragmentation 
mass is a strong function of metallicity (see Fig. [8| . At 
low metallicitics (Z < 10~ 4 Z Q ) one large fragment is 
produced. Conversely, high metallicities (Z > 1O _4 Z ) 
produce ~ 10 3 — 10 4 solar mass fragments, lending sup- 
port to metallicity being the driver of the Pop III-II 
star formation transition protogalactic accretion flows. 
This behavior can be understood by comparing the cool- 
ing time from metals and a c hara cteristic fragmentation 
time, as we explain in Section |4.3| below. We note that in 
runs with PISN abundances, a given metallicity will pro- 
duce fragments nearly one order of magnitude smaller. 
Yields from PISNe are thus more effective at reducing 
the fragmentation mass scale. 

4.2.1. The Impact of Dust 
Th e presence of dust, whic h may accompany metals 



(e.g., Cazaux & Tielens 20021, could have a significant 



impact on the thermal evolution of the gas in that it ra- 
diates thermally and also catalyzes molecular hydrogen 
formation. Dust grains should be produced in primordial 
supernovae dTodini fc Ferrara||2001| |Nozawa et al~1|2003 ' 



Schneider et al. 2004), but uncertainty still exists as to 
ho w much dust ultimately enriches the ISM (Nozawa et 
aLl |2007| |Cherchneff fc Dwek|l2010[ ) . We do not include 
dust thermal emission or gram surface H2 and HD for- 
mation in our chemical model. Dust thermal emission 
is always a sub-dominant coolant at met allicities up to 



10 Zp, and densities below 10 cm d (Omukai et al. 



2005 1 . To assess the catalytic impact of dust grams m our 
model, we compare the H 2 formation rate on dust grains 
with the H2 formation rate through the BP channel and 
the LW H2 photo-disso ciation rate. The format ion rate 
of H2 on dust grains is ( Cazaux fc Spaans||2004 ) 



i? d = 3.025 x 10- 17 e H2 n H n(H) 5 H (T) / a (T) 
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Fig. 8. — Characteristic fragmentation mass as function of metal- 
licity for three different metallicity abundance patterns. The 
masses represent the Bonnor-Ebert mass at the point of fragmen- 
tation, when i soun d = tg. These represent the physical scenario 
m case A (J21 = 10 4 ) but with varying levels of metallicity. The 
two curves labeled PISN (overlaid) are for a 150 and 250 Mq PISN 
mctallicity yield pattern. The critical metallicity (the metallicity 
where the fragmentation mass scale drops appreciably) occurs at 
slighly lower metallicity for the two runs performed with PISN 
abundance patterns, due to the higher proportion of oxygen and 
silicon produced in PISNe as compared to the solar abundance pat- 
tern, the shift is small however. The upper mass scale of ~ 10 5 Mq 
reflects the initial conditions of the accretion flow while the lower 
mass limit of ~ 3 Mq i s set by the CMB temperature floor. As 
discussed in Section |4.3| the critical metallicity can be understood 
by considering the three relevant timescales, t coo i, t soun d, an d *ffi 
it is dependant on the gas being able to cool appreciably before the 
fragmentation criteria is reached (determined by t aoun d and iff ) . 



where eH? i s the H2 formation efficiency (Cazaux fc 
Spaans|2009 ), S'h(T) is the temperature depe ndent stick- 



ing factor o f hydrogen onto dust grains (Cazaux fc 
Spaans|[2004l , Tin is the number density of neutral fry 



drogen, n(rij is the number density of hydrogen nuclei, 
and f&(T) is a temperature dependent factor account- 
ing for t hermal evaporation of hydr ogen from the grain 
surface (Tielens fc Hollenbach 1985). We assume a typ- 
ical dust grain has a mass of 1.2 x 10 -14 g, such that 
rid = 1-3 x 10 _10 ^rfnH- Equation (10) is normalized for 
a standard Galactic dust-to-gas mass ratio of £d = 0.01 
at around solar metallicity. We can scale this for lower 
metallicity by writi ng (Z) = (Z/Zq) R^Zq) where 
i?d(Z©) is equation (10) evaluated for £d = 0.01. 

In Figure [9] we compare these three volumetric chemi- 
cal rates — the LW H2 destruction rate, the H2 formation 
rate through the B~ channel, and the H2 formation rate 
on dust grains for two different values of Jji- The grain 
H2 formation rate is computed with equation ( 10 ) using 
the isobaric density evolution scheme while the other two 
rates are extracted from our chemical model. In the run 
with J 2 i = 10 4 , LW photon destruction of H 2 dominates 
over dust grain H2 formation at temperatures T > 250 K. 
Below this temperature, H2 formation on dust grains is 
extremely efficient. We also repeat this run for J21 = 10 
and find similar behavior. We chose these two values 
of J21 as they bracket the critical J 2 i for H 2 formation 
through the B~ channel (see Section 4.2). We do not 



(10) consider dust opacity and its effect on the LW photon 
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Fig. 9. — The volumetric chemical rates for the LW H2 destruc- 
tion rate (labeled LW), the formation rate of H2 through the H - 
channel (H — ), and the grain H2 formation channel rate (Dust) for 
two different values of the intensity in the LW bands (see Section 
|3.3| ; J21 = 10 4 (top) and J21 = 10 (bottom). The dust forma- 
tion rate drops off rapidly at high temperatures due to the steep 
temperature dependences of the H sticking coefficient Su (T) and 
the H evaporation coefficient / a (T). In the high LW intensity run, 
LW destruction initially dominates the chemistry of H2 while dust 
formation is negligible. At 600 K, H2 is in equilibrium with the 
radiation field and at 300 K dust H2 formation becomes the pri- 
mary formation channel. Gas which reaches roughly T ~ 250 K 
will experience rapid H2/HD formation, allowing the gas to cool 
rapidly to Tomb- 



flux, as understanding the interplay between photo-disso- 
ciation and dust requir es complicated multidimen sional 
modeling of PDRs (e.g., |Krumholz et al.|2008||2lj09"} . We 
also do not consid er H2 formatio n heat ing in this work, 
but as is shown in Omukai et al. (2005), at metallicities 
up to 10~ 2 Zq and densities below 10 b cm _ 3, this heat- 
ing mechanism is subdominant compared to H2, metal 



line cooling, and adiabatic heating. 

Examining Figure [8] and using equation ([£}, T ~ 250 K 
corresponds to a BET mass of ~ 230 M which is the 
initial fragmentation mass scale of a gas with metallic- 
ity Z ~ 10~ 3 Z Q . Molecular hydrogen, however, is not 
able to cool gas below T ~ 200 K due to its lack of a 
permanent dipole moment, thus further cooling depends 



on the formation rate of HD and other molecules (e.g., 
CO) which can cool the gas to even lower temperatures — 
these molecules will most likely also be formed on dust 
grains in significant number around this temperature as 
well. We can expect then, because of dust's impact on 
molecule formation, that the fragmentation mass scale 
will be sharply lowered at around 10~ 3 Z Q . At higher 
mctallicity ~ 10~ 2 Zq, the formation rate of H2 on dust 
grains will increase. However, gas with a metallicity of 
10~ 2 Zq will reach Tqmb through metal cooling alone 
(see Figure [8}; the addition of molecular cooling will only 
decrease the cooling time. Thus, while the formation 
rate of H2 on dust grains will become increasingly larger 
at higher metallicity, it will not affect the end product, 
namely the initial fragmentation mass scale, of the gas in 
our model. Increased molecule formation brought about 
by dust may also allow thermal instabilities to develop 
more easily — see Section |2~4j We conclude that molecule 
formation on dust grains will alter the behavior of the 
gas below some critical temperature (T ~ 250 K) mak- 
ing the transition from high mass fragments to low mass 
fragments occur more abruptly, but while retaining the 
general picture which we present here. 

These find ings are complementary to those of [Omukai 
et al. ( 2008 ) who also use a one-zone model to study the 
thermal evolution of gas in a halo with T vir - 10 4 K. 
They find, with the inclusion of dust, that the gas will 
fragment vigorously with a metallicity as low as Z ~ 
10 -6 Zq. This episode of dust- facilitated secondary frag- 
mentation occurs at densities uh ~ 10 10 cm -3 , higher 



than the densities, 



10 



5-6 



cm , at which we first 



encounter fragmentation. For the densities which our 
models overla p there is exce l lent a greement — see Figures 
5a and 5b in Omukai et al. ( 2008 ) . 



4.2.2. Initial Density 

These results depend sensitively on the choice of initial 
post-shock density, due to the dependence of th e cooling 
rate and the BE mass on density. Motivated by IGreif et 

* I J ~L . . _ L_5 — j 



al. ( 2008 ) , we chose a pre-shock density of n = Hrcm 



as this is the density right before the cold stream gas 
experiences a shock near the center of the galaxy. If 
the pre-shock number density is smaller, 10 2 cm -3 , the 
initial BE mass is ~ 5 x 10 5 M@. The critical metallicity 
for the gas to reach Tqmb is also correspondingly shifted 
to higher metallicity roughly 10 -2 5 Zq, and even at this 
metallicity, the gas fragments with Mbe ~ 10 M . For 
a pre-shock density of 10cm~ 3 , the initial BE mass is 
~ 2 x 10 6 Mq and the gas temperature is not able to 
reach Tcmb at any reasonable metallicity. 

4.2.3. Hot Accretion Mode 

We proceed to comment on the implications of an even 
lower pre-shock density encountered in quasi-spherical 
baryonic accretion from the IGM onto the halo . As has 
been shown in simulations (iWise & Abel 20071 IGreif et 



al.|2008 ), accompanying the cold accretion mode is a hot 
mode where gas is accreted at the virial radius directly 
from the IGM residing between the (cold) filaments of 
the cosmic web. The fate of gas accreted in this mode 
will differ substantially from the case of cold accretion, 
mainly in the location of the shock and the subsequent 
cooling rate. As gas is accreted directly from the IGM, 
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its post-shock density is ~ 0.4 cm" 3 ( [Wise fe Abel|2007| ), 
much lower than the initial density adopted tor this study 
(10 3 cm -3 ). Initially, the gas is not self gravitating as 
indicated by its large BE mass. None of the cooling pro- 
cesses discussed here will allow the gas to cool appre- 
ciably below T ~ 7000 K. At this point, t cooi ~ 10 9 yr 
while t g ~ 3 X 10 7 yr; this suggests a stable virial shock 
may develop, even with a high degree of metallicity in 
the gas. It seems likely that only the gas which enters 
the galaxy through a dense, cold accretion filament is 
able to immediately fragment with a small characteristic 
fragmentation mass. 

4.3. Analytical Toy Model for Critical Metallicity 

As shown in Section |4.2[ runs with a large LW radia- 
tion background ( J21 > 10 ) exhibit a strong dependance 
on metallicity of the final fragment mass scale. Here, 
we discuss the physical basis behind this critical metal- 
licity, utilizing a simple toy model. As demonstrated, 
our critical metallicity is a result of cosmological initial 
conditions via the virial velocity and accretion density, 
cooling physics, and the Jeans analysis for gravitational 
fragmentation. 

Let ifr a g be the time it takes the expanding shocked 
region to reach the fragmentation criterion: i S ound > iff; 
where t soun( j ~ L/c s , the ratio of the length of the 
shocked region, given by equation Q, and the sound 
crossing time. In the absence of metals and molecular 
cooling, atomic cooling quickly brings the temperature 
from its post-shock value of T ~ 10 4 K to the asymptotic 
T ~ 6,500 K, resulting in Lz=a{t) ~ | v^t, where the 
factor of 1/5 is from equation ([8| considering adiabatic 
and isobaric density enhancement at T — 6,500K. In 
metal- and molecule-free gas fragmentation consequently 
sets in at if ra g,z=o ~ 10 6 yr. Interestingly, due to the 
competing temperature dependence of the free-fall and 
sound crossing time, if rag is roughly constant for any 
value of metallicity. To determine the level of metals 
needed to allow the gas to fragment at a small mass scale, 
we evaluate tf rag ~ t coo i, where £ coo i is the cooling time 
due to metal transitions. Equating the timescales yields 



tc 



t{Z) 



-A 10 AE 1Q ~ 10 6 yrs, 



(11) 



where A w is the spontaneous emission coefficient of a 
transition, AE\q is the corresponding energy difference 
between states, and X\ = ni/n to t is the number density 
fraction of the metal coolant. For the solar abundance 
pattern, O is the dominant coolant over a wide range of 
temperatures (see Fig. |2J). For the two lower states of O, 



A w = 8.9 x 10~ 5 s" 1 and AE W = 3.14 x 10" 
critical O abundance therefore is 



knT 



10 6 yrs AmA^i 



10" 



'ergs. The 



(12) 



where T — 6, 500 K. This fraction corresponds to an oxy- 
gen mass fraction Zq ~ 10~ 48 Zo,©; agreeing with the 
critical metallicity demonstrated by Figure [S] The same 
analysis could be performed with the cooling rate of sil- 
icon (as it is the most effective coolant for runs with a 
PISN abundance pattern) — the conclusion would remain 
the same. Previous studies on the role playe d by metal- 
licity in the star formation mode transition ( |Bromm et| 



aT] |2001b[ pmukai et aT1|2005l [Santoro fc Shull| 20061 
have defined the critical metallicity as the metal frac- 
tion that alters the characteristic mass of stars formed 
from the standard Pop III initial conditions where cool- 
ing by H 2 and metals competes with the compressional 
heating due to self-gravity. Our physical scenario is very 
different; instead of having a self-gravitating, adiabatic 
collapse in 10 6 Mq DM minihalos, we consider an isobari- 
cally cooling gas that was shocked to 10 4 K in 10 8 Mq ha- 
los. Direct comparison between critical metallicity values 
obtained in these two contexts is not entirely justified. 
Nonetheless, we do find a critic al metallicity compara ble 
to some of the minihalo results ( |Bromm et al.|2001b[ ). It 
is not clear why a hot, isobarically cooling shock should 
behave so similarly to the adiabatic collapse experienced 
in a minihalo. 

5. DISCUSSION AND CONCLUSIONS 

In this work, we have examined the fate of metal- 
enriched, cold, dense, and dust-free accretion streams 



into DM halos of mass 10 s M m at redshift 



10; we 



regard these objects as the first galaxies. From simula- 
tions, we know these cold streams penetrate very close 
to the center of the protogalaxy, where they shock and 
fragment. Our goal was to determine the fragmentation 
properties of these cold streams and to identify the key 
physical parameters governing this first step in the high 
redshift star formation process. We find that metallic- 
ity (due to metals being an effective low temperature 
coolant) and LW background strength (as it is capable 
of destroying H2 and HD) are the dominant factors in de- 
termining the final fate of the post-shock region. Using 
a fit to the baryonic growth rate of halos, and assuming, 
as recent studies suggest, that these galaxies are fed en- 
tirely by cold streams and that their star formation rate 
closely matches the mass accretion rate, we can estimate 
the number of fragments produced by these flows. We ei- 
ther see one large (~ 10 4 M o ) fragment being produced, 
or roughly 10 4 solar-mass fragments, depending on the 
initial conditions of the flow. The CMB temperature 
floor is also a key factor in determining the mass of frag- 
ments, as it represents a strict temperature lower limit 
of radiatively cooling gas. Roughly speaking, this sce- 
nario occuring at higher redshifts would still demonstrate 
similar behavior in terms of a critical metallicity, but 
would produce fragments of slightly higher mass. Over- 
all, this effect is of little importance as even low-redshift 
star forming molecular clouds do not generally cool be- 
low 10 K, wheres as the CMB temperature at z = 10 is 
~ 30 K. Any observational effects from this temperature 
floor wold be difficult to extract from our simulations. 



Bailin et al. (2010) studied CMB regulated star forma- 
tion m some detail (at lower redshifts than we examine) 
and find that observations currently cannot confirm or 
refute CMB regulation on star formation. 

We did not include dust, even though it can have a 
major impact in that it allows very efficent H2 forma- 
tion on grain surfaces, though we do argue its presence 
will not significantly alter our findings. The presence of 
dust and its effect on cooling has been discussed as a 
possible catalyst for th e sta r formation mode transition - 
Schneider etaE] ( |2006[ ) and jSchneider fc Omukai| ( |2009 1 
found that dust induced fragmentation can decrease the 
Pop II fragmentation mass scale by about four orders of 
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magnitude over that accessible wi th metal fine structure 
lines only. Jappsen et al. ( 2009b I suggested that metal- 
licity does not attect the fragmentation of gas with hot 
initial conditions, similar to the environment probed in 
this study. We attribute the difference in our conclusions 
to our examination of a slightly different physical regime 
and to our inclusion of a soft UV b ackground that sti - 
fles H 2 cooling (as supported, e.g., by Shang et al.|2009 ). 
Taking guidance from galaxy formation simulations, we 
follow isobarically evolving gas in an expanding shocked 
slab to study its thermal evolution and fragmentation. 
Our model is designed to reflect the filament-dominated 
inflow topolog y of cold accretion flows, whereas [Jappsen] 



et al. (2009b) examine the collapse into a spherically- 



symmetric, centrally concentrated halo potential. In our 
isobarically evolving case, the tight coupling between 
density and temperature will increase the importance of 
the role of the cooling rate in determining the final out- 
come of fragmentation. This aside, if we include molecule 
formation and molecular coolin g (H2 and HP) by sett ing 
J 21 = 10, our results agree with Jappsen et al. (2009a I in 
that metal line cooling does not attect the fragmentation 
properties of the evolving gas. This emphasizes the im- 
portance of a LW soft UV background in modeling the 
formation of Pop II stars in the first galaxies. 

The first galaxies are suspected to be the first struc- 
tures where turbulence begins to play a major role, 
driven by high Mach-nu mber cold streams penetrating 



to th e protogalactic core (Greif et al.|2008 Wise & Abel 



20071. It is the interplay ot turbulence and sett grav- 



lty which is suspected to be responsible for the shape 
of the present-day IMF, especially the high mass power 
law extension (Padoan & Nordlund 2002; Mac Low & 
Klessen 2004; McKee & Ostriker 2007). Thus, it is likely 
that star formation inside the first galaxies will be char- 
acterized by a stellar mass distribution that for the first 
time in the history of the universe resembles already the 
present-day IMF. Specifically, we can expect turbulence 
to introduce density perturbations on the order of M 2 , 
suggesting within our above model that fragments as low 



as ~ 0.3 M & might form. Such stars, below the mass 
limit of 0.8 M Q required to survive for the age of the uni- 
verse, would still exi st today, and could be probed by 



Stellar Archaeology (Frebel et al. 2009 and references 



therein). To predict the detailed shape of the Pop. II 
IMF inside the first galaxies, the source and impact of 
the turbulence needs to be studied with realistic numer- 
ical simulations. 

Here, we have constrained some of the key physical 
processes that govern star formation in high-rcdshift pro- 
togalaxies within an idealized analytical framework that 
allowed us to explore the relevant parameter space. This 
work needs to be followed up by detailed simulations, 
including a number of ingredients ignored here, such 
as hig h optical depth effe cts, H and He ionizing radi- 

feedback effects 
an d possibly cosmic rays 



2002}, turbulence 
stellar outflows 



etc. 



ation (Ricotti et al 
(AGN, 

from the first SNe ( Jasche et al. 2007| |Stacy fc Bromm] 
2007|). It needs to be emphasized that the estimates pre- 



sented in this paper are derived with the Bonnor-Ebert 
mass and a simplified treatment of the post-shock dy- 
namics. To truly probe the fragmentation mass scale in 
the first galaxies one must use multidimensional numer- 
ical simulations. The goal of understanding the com- 
plex first stages of galaxy formation seems now within 
reach, using state-of-the-art three dimensional radiation- 
hydrodynamical simulations, together with realistic cos- 
mological initial conditions. Such calculations are crucial 
to make useful predictions for the James Webb Space 
Telescope (JWST). 
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TABLE 3 Reaction Rates for our Chemical Network. 



Reaction 




Rate Coefficient (cm 3 s 1 ) 


j Notes Reference 


H + c~ -> H~ +7 


fcl = 
+ 


dcx[- 17.845 + 0.762 logT + 0.1523(logT) 2 

U.UOZ l^lOg 1 1 J 

dcx[-16.420 + 0.1998(logT) 2 
5.447 x 10" 3 (logT) 4 
4.0415 x 10- 5 (logT) 6 ] 


1 ^ 0, UUU rv 
T > 6.000K 


H- +H -> H 2 + e" 


fc 2 = 


1.5 x 10~ 9 

4.0 x 10~ 9 T~°' 17 


T < 300 K 
1 > 300 K 


H + H+ H+ +7 


fc 3 = 
+ 


dcx[-19.38 - 1.523 logT 
1.118(logi ) — 0.12o9(logi )°J 




H + H+ -)■ H 2 + H+ 


&4 = 


6.4 x lO" 10 




H" +H+ H + H 


fc 5 = 


5.7 x io- 6 T-°- 5 + 6.3 x 10- 8 
9.2 x 10- n T°- 5 + 4.4 x 10~ 13 T 




H+ + e- -> H + H 


fc 6 = 


1.0 x 10" 8 

1.32 x i(r 6 T- - 76 


T < 617K 
T > 617K 


H 2 + H+ -> H+ + H 


fc 7 = 


[-3.3232183 x 10~ 7 
+ 3.3735382 x 10" 7 lnT 

- 1.4491368 x 10~ 7 (lnT) 2 
+ 3.4172805 x 10~ s (lnT) 3 
-4.7813720 x 10' 9 (lnT) 4 
+ 3.9731542 x 10~ 10 (lnT) 5 

- 1.8171411 x lO-^lnT) 6 
+ 3.5311932 x 10" 13 (lnT) 7 l 
xexp(- 212 T 37 - 15 ) 




H 2 + e~ -» H + H + e" 


fc 8 = 


3.73 x 10- 9 T 01121 exp(- 99430 ) 




H 2 +H^H + H + H 


^9,1 
™cr,H 


= 6.67 x 10- 12 T°- 5 cxp[ (1.0 + 63590 )] 
= 3.52 x 10- 9 exp( 43900 ) 

= dex[4.0 - 0.416 log (^g) - 0.327 log 2 (^j 




H 2 + H 2 -> H 2 + H + H 


felO.l 
felO.h 
n cr,H 2 


5.996xl0- 30 T 4 ' 1881 / 54657.4-, 
~~ (1. 0+6. 761xlO- 6 T)5 6881 V T > 

= 1.3x 10- 9 cxp( 53300 ) 

= dex [4.845 1.3 log ( J m) ) + 1.62 log 2 ( ^ )] 




H + c~ — > H+ + c~ + c~ 


fell = 


= cxp[-3. 271396786 x 10 
+ 1.35365560 x 101nT c 

- 5.73932875 x 10° (lnT e ) 2 
+ 1.56315498 x 10° (lnT e ) 3 

- Z.o ( f UodUU X 1U (in 1 e ) 
+ 3.48255977 x 10~ 2 (lnT e ) 5 

- 2.63197617 x 10- 3 (lnT e ) 6 
+ 1.11954395 x 10" 4 (lnTe) 7 
-2.03914985 x 10- 6 (lnT e ) 8 ] 




D + c~ -> D+ + e~ + c~ 


fel2 


= fen 
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TABLE 3 — continued from previous page 





Reaction 




Rate Coefficient (cm 3 s i) 


1 Notes 


Reference 


H+ + e" 


-> H + 7 


fel3 


= 2.753 X 10-14(315614)1.503 
x[L0 + (604625)0.4701-1.923 






D+ + e" 


-)■ D + 7 


fel4 


= fcl3 







k 15 = cxp[-l. 801849334 x 10 

+ 2.36085220 x 10°lnT e 
-2.82744300 x 10 _1 (lnT e ) 2 
+ 1.62331664 x 10~ 2 (lnT e ) 3 
H" + c~ H + e- + e- - 3.36501203 x 10~ 2 (lnT e ) 4 

+ 1.17832978 x 10- 2 (lnT e ) 5 

- 1.65619470 x 10- 3 (lnT e ) 6 
+ 1.06827520 x 10- 4 (lnT e ) 7 

- 2.63128581 x 10- 6 (lnT e ) 8 ] 



H-+H^H + H + e" 




= 2.5634 x 10- 9 T c i' 7S i 86 
= exp [-2. 0372609 X 10 

+ 1.13944933 x 10° (lnT e ) 

- 1.4210135 x 10-1 (lnTe) 2 
+ 8.4644554 x lO -3 (lnT e ) 3 

- 1.4327641 x 10- 3 (lnT e ) 4 
+ 2.0122503 x 10- 4 (lnT e ) 5 
+ 8.6639632 x 10- 5 (lnT e ) 6 

- 2.5850097 x 10 (lni e ) 
+ 2.4555012 x 10- 6 (lnT e ) 8 
-8.0683825 x 10- 8 (lnT c ) 9 ] 


T c < 0.1 cV 
T c > 0.1 cV 


H~ +H+ -»■ H+ + e" 


fcl7 


= 6.9 x io- 9 T-°- 35 
= 9.6 x 10-7T-0-90 


T < 8000 K 
T > 8000 K 


H + D+ D + H+ 


fcis 


= 2.06x 10-i°T - 396 exp(-f ) 
+ 2.03 x i - 9 T-°- 332 




H+ + D -> H + D+ 


fel9 


= 2.0x 10-i°T°' 402 cxp(- 37 -i) 
- 3.31 x iO"i 7 Ti- 48 
= 3.44 x io-i°T°- 35 


T < 2 x 10 5 K 
T > 2 x 10 5 K 


H 2 + D+ -> HD + H+ 


fc 2 () 


= [0.417 + 0.846 logT - 0.137 (logT) 2 ] x 10" 9 




HD + H+ -»■ H 2 + D+ 


fc 2 l 


= 1.1 x 10- 9 cxp(-^) 




D + H 2 ->■ HD + H 


&22 


= 1.69 x 10-1° exp (- 4 f°) 

= 1.69 x 10-1° exp ( 4 ff 8 ° + iiff°°) 


T < 200 K 
T > 200 K 


HD + H -> D + H 2 


&23 


= 5.25x 10-n exp (- 4 f°) 

= 5.25 x 10-n exp ( ifo + iTpoO) 


T < 200 K 
T > 200 K 


He + e- -> He+ + e~ + e~ 


&24 


= exp [-4.409864886 x 10 
+ 2.391596563 x 101nT e 

- 1.07532302 x 10(lnT o ) 2 
+ 3.05803875 x 10° (lnT e ) 3 
-5.6851189 x 10-i (lnT e ) 4 
+ 6.79539123 x lO" 2 (lnT e ) 5 

- 5.0090561 x 10- 3 (lnTe) 6 
+ 2.06723616 x 10- 4 (lnT e ) 7 
-3.64916141 x 10- 6 (lnT e ) 8 ] 





Continued on next page 



FRAGMENTATION IN THE FIRST GALAXIES 17 

TABLE 3 — continued from previous page 



Reaction 




Jrtate Loemcient (cm° s ) 


1 JNotes Reference 


He+ + e" -> He + 7 


fc25,rr 


= io- n T-°- 5 [11.19 - 1.676 log T 

— U.zooz (log _Z J + U.U44oo(lOgi ) \ 

= 1.9x 10- 3 T- 15 exp( 473 T 421 ) 
x [1.0 + 0.3 cxp( 9 « 84 )] 




He+ + H -> He + H+ 


&26 


= 1.25x 10- 15 ( 3^5 ) ' 25 




lie + ri T — ► lie 1 + ri 


&27 


= 1.26x 10- 9 T-°- 75 cx P ( 127 T 500 ) 
= 4.0 x i(r 37 T 4 - 74 


T < 10 4 K 
T > 10 4 K 


He+ + D -> He + D+ 


k 2S 


= fc 26 




He + D+ -> Hc+ + D 


&29 


= k 27 




C+ +c" -> C + 7 




= 4.67xl0- 12 (3^)~°' 6 

= 1.23xl0- 1 7 (^) 2 ' 49 exp( 2 l^) 

= 9.62xl0-«(^)- 1 ' 37 exp(-i™ ) 


T < 7950 K 

7950 K < T < 21140 K 
T > 21140 K 


Si+ + c~ -> Si + 7 


k;u 


= 7.BxlO-«(^)-° o " ffl 
= 4.86xl0- 12 (3^) ' 32 

/ n 0.818 

= 9.08xl0- 14 (^) 


T < 2000 K 

2,000if < T < 10,000 K 

i > 1U, UUU xV 


0+ + e~ O + 7 


k-i2 


= 1.30 x io- 1() T-°- 64 

= 1.41 x io- 10 T-°- 66 + 7.4 x 10- 4 T- 15 

xcxp( 17 « K, °) [1.0 + 0.062 cxp( 145 T 000 )] 


T < 400 K 
T > 400 K 


C + c~ -> C+ + e" + e~ 


k:i3 


= 6.85 x 10- 8 (0.193 + u)- 1 u°- 25 e- u 


u = 11.26/T 


Si + c- -> Si+ +c" +e" 


k 3 4 


= 1.88 x 10~ 7 (1.0 + m°' 5 )(0.376 + u)- 1 u°- 25 e~ u 


u = 8.2/T c 


+ c~ -> 0+ +e- +e~ 


k-i5 


= 3.59 x 10- s (0.073 + u)- 1 u°- 34 e~ u 


u = 13.6/r e 


0++H^O + H+ 


&36 


= 4.99 x io-Ht - 408 + 7.54 x io- 10 T-°- 458 




O + H+ ->■ 0+ + H 


^37 


= [1.08 x io- n T - 517 

+4.00 x 10- 10 T - (,0669 ]cxp(- 2 |?) 




O + Hc+ -> 0+ + He 


k-S8 


-1991"10~ 15 f T ~)°' 3794 crDf T ) 
— 4.yyi x iu (.loooj)-' p[ - 1121000 > 

+2.780 x 10- 1B ( 10 ^ 00 )- - 2163 ex P ( 81 f 800 ) 




C + H+ -> C+ + H 


&39 


= 3.9 x lO-iST 213 




C + H+ -> C+ + H 


fc 40 


= 6.08xl0-i 4 (^)i- 9 6ex P (-l 7 f»«) 




C + Hc+ ->■ C+ + He 


k 41 


= 8.58 x io- 17 T - 757 
= 3.25x io-^T 968 
= 2.77 x io- 19 r 1597 


T < 200 K 

200if < T < 2000 K 

T > 2000 K 
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TABLE 3 — continued from previous page 






Reaction 




Rate Coefficient (cm 3 s — 


Notes 


Reference 


Si + H+ -> Si+ + H 




= 5.88x i - 13 T - 848 
= 1.45 x 10~ 13 T 


T < 10 4 K 
T > 10 4 K 




Si + Hc+ -»■ Si+ + He 


k43 


= 3.3 x 10" 9 






C+ + Si -)■ C + Si+ 


fc44 


= 2.1 x 10" 9 








^75,1 


= dcx [-27.029 + 3.801 logT - i^ 1 )] 






H 2 + He -> H + H + He 


fc 75,h 
1cr,Hc 


= dcx[ 2.729 1.75 logT ( 23474 )] 
= dex [5.0792 - 6.247 x 10~ 5 (T - 2000)] 






H 2 + He+ -> He + H+ 


fcl3S 


= 7.2 x 10~ 15 






H 2 + He+ -> He + H + H+ 


^139 


= 3.7 x 10- 14 cxp(-|l) 






H + H + H^H 2 + H 


fc 2 06 


= 1.32x 10" 32 (^)-°- 38 
= 1.32 x IO- 32 (ggj)" 10 


T < 300 K 
J > 300 K 




H + H + H 2 -> H 2 + H 2 


&207 


= 2.8 x 10~ 31 T~ ' 6 






H + H + He -> H 2 + He 


fc 2 os 


= 6.9 x icr 32 T-°- 4 






Fe + e~ -> Fc + + e~ + e" 


fc 222 


= 2.52 X 10~ 7 (0.701 + u)" 10 0)°' 25 exp(-u) 


u = 7.9/T 


1 


Fe+ + e- Fe + 7 


fc 22 3 


= 1.42x10 (3^1) 

+ 1.6xl0- 3 T- ls exp( 60 * 1()4 ) 
[l.O + 0.46 cxp( 8 - 97 ^ 104 )] 




2 


He+ + e" -> He++ + e~ + e" 


^300 


= 5.68 x 10- 12 cxp( 631 T 515 ) 
x (T~ °- 6 4- 3 16228 x in -3 ") -1 




3 


He++ + c~ -> Hc+ + 7 


fool 


= 1.34 x 10- 9 T-°- 7 [1.0 + (^s) - 7 ]- 1 - 




3 



